Coefficient functions and open charm production in deep inelastic 

scattering 

A.V. Kisselev* 
Instituto de Fisica Gleb Wataghin, UNICAMP 
13083-970, Campinas, SP, Brasil 
and 

Institute for High Energy Physics, 142284 Protvino, Russia^ 





Abstract 

It is shown that the problem of double counting in open charm production 
in DIS can be solved by using the expression for DIS coefficient functions in 
terms of 2PI diagrams. 

PACS number(s): 13.60.Hb, 13.85.Ni, 14.65.Dw, 12.38.-t 



*Email address: kisselev@mx.ihep.su 
^ Permanent address 



1 



Open charm production in deep inelastic scattering (DIS) is a subject of interest from 
both experimental and theoretical point of view. Recent data from HI ]I[ and ZEUS |2j 
collaborations have shown that charm quark contribution is an important component of DIS 
structure functions. 

On of the most predictive methods of calculating open charm contribution to the struc- 
ture functions is fixed order perturbative QCD (see ||). Much efforts have been also done 
in order to formulate variable flavor number scheme (see 0-0]). 

According to the factorization theorem P,p|,|lO|,p!l| , the contribution of charm, F 2 C , can be 
represented as follows^: 



-F^Q\x)= 

X Jx Z 



C\Q\ /i 2 , z)q c ^\ -) + C 9 (Q 2 , /i 2 , z)g^ 2 

z z' 



(1) 



Here C a are the process-dependent coefficient functions (a = c, g) , q c and g being charm 
and gluon densities of the incoming hadron and /i is a factorization scale. The gluon coef- 
ficient function, C 9 , includes, in particular, photon-gluon fusion (PGF) contribution which 
dominates at low scales, Q 2 < m 2 , where m c is a mass of charm quark. At high Q 2 , as was 
noted in Refs. [f|||, part of the PGF cross section is generated by the evolution of the charm 
contribution (the first term in Eq. ([!])). 

To avoid double counting, one has to subtract this contribution. In Ref. [||] this subtrac- 
tion has been taked into account in lowest order in a s (in what follows a symbol <g> means 
a convolution in variable z): 

F 2 = Cf } ® q c - Cf > ® P cg <g> g + C« <g> g + 0(a 2 ) . (2) 

If fj, ~ m c , we have an approximate cancellation of the first two terms in Eq. (0) and we 
arrive at the dominance of PGF mechanism. On the other hand, if fi 3> m c , the last two 
terms in Eq. (Q) are almost cancell and we reproduce the QCD parton model in the leading 
order. 



1 Let us note that at very small values of x one has to consider more general k± -factorization [12 
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In fact, in Eq. (|2]) one should add another order a s contribution of the form 

CP ®q c ~ <g> P dc ® q c (3) 

which is numerically less important. 

Let us note that the Wilson type coefficients in ([I]) are different from the coefficient func- 
tions C c and C g in (0) and (|3]). While the formers have no infrared or collinear singularities, 
the same is not true for the latters. 

In Ref. H the modification of the gluon coefficient function has been proposed: 

Cp - C* = C ™ F ~ AC *- ( 4 ) 

By using the well-known PGF expression in the first order in a s (see, for instance, |JL3[|), 
ACW was calculated in @. At Q 2 /m 2 c > 1 it looks like 



Ztt 



o 2 

P cg {z)\n\ + z{l - z) 



(5) 



The corresponding order a 2 expression was also found (with the account of only large loga- 
rithmic terms) ||: 

ACf = C?> ® [Kin Q 2 ) 2 (P£> ® P 9 ? + Pj? ® F<») + a>g 2 pw; 

+ C« ® a s In Q 2 P£> + a s C'( 1) ® a s In Q 2 P$. (6) 

Both procedures, (0) and (Q), treat the problem of double counting by hand and they 
cannot be easily generalized to higher orders in a s . Let us note that expression @ was 
obtained with the account of only leading logarithmic contributions in each subtracted term, 
and the factorization scale was assumed to be large, fi 2 ~ Q 2 . 

To overcome these shortcomings, let us start from a definition of coefficient functions 
in terms of two particle irreducible (2PI) amplitudes. In the following we assume that a 
corresponding object is 2PI in the direction of the iteration (that is in the t-channel). We 
will work in an axial gauges (n^A^ = 0, being a gauge vector) and follow a scheme 
developed in Ref. II. 



Let D being a matrix of formal parton distributions (below we define physical parton 
distribution functions D (0)). It obeys a matrix equation 

D= (I + V + V 2 + ■■■)£>= (I-V)- 1 !), (7) 

where D is 2PI part of D, 2PI kernel V defines the evolution of D and J is a unit matrix. 
To simplify notations, we dropped a sum both in parton types and spins. Integration in 
internal momenta are also assumed in (fj). 

It is known that in the axial gauge 2PI aplitudes have no singularities associated with 
a propagation of intermediate physical states f8|. For instance, V(k n , k n -\) in © is finite 
at k\_ x — > 0, but it includes legs, corresponding to momentum k n , and has a pole in k\ as 
A£-0. 

Then DIS structure function has the form 

F = F + A(I -V)- 1 !) (8) 

with F being the 2PI part of F, while A is the 2PI part of the virtual photon-parton 
amplitude. 

Following 0, let us introduce a projector operator onto physical helicity states (d^ u is a 
tensor part of the gluon propagator) 

(k) a 13 ] [ {—^—)ai3i for a quark line 
i Akn (9) 

-gP " ] [ (-g^), for a gluon line 
and another operator which projects onto small virtualities of parton linesf] 

p, = e^-\k 2 \). (io) 



2 Another possibility is to use an operator P £ which extracts poles in e in MS renormalization 
scheme (see, for instance, fl4jl). We prefer to use the operator in order to have clear physical 
meaning for the scale \i (see below). 
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As we will see, our result ( j22|) depend on a product of Ph and P^ (p3|). But at intermediate 
steps it is convenient to consider these operators separately. 
By using matrix identities 



A-l 



(/ - v)- 1 = [/-(/- P h )V)]-\l - P h V 
(I - PhVT 1 = (I ~ V)- l [I + P^P h V(I - V)' 1 



with 



V' = V[I -(I- P h )V]-\ V=(I- P,)P h V, 



we can rewrite Eq. @ in the form 

F = A(I-V)- 1 
+ A(I-V)- 1 

where 



k i =fc_=0 



k i =fc_=0 



£> + Ai(J - AV)" 1 /) 
Ay(J - Ay)" 1 /) + F, 



i = i[/-(j-p ft )y]- 1 



and 



AA = A-A{q,k) 
AV = V-V(r,k) 



k I =fc_=0 



fc I =fc_=0 



Here fc_ = (k 2 + k\)/kn. 

In Eq. (|13D we defined the quantity 



D = [I + P^P h V(I-V)- 1 ]D. 



(11) 



(12) 



(13) 



(14) 



(15) 



(16) 



The quantity AA(I — AV)" 1 in the second term in Eq. ( |13"D has no contributions of the 
type Q!g(lnQ 2 ) . As far as one is interested in logarithms of Q 2 , it is possible to omit this 
therm in the coefficiet function. The third term in Eq. (|T3"D, AV(I — AV)~ 1 D, is of the oder 
of A 2 //i 2 , where A is a tipical hadron scale. The non-perturbative scale A arises as a result 



of integration of the kernel V with D (which describes initial parton distributions inside the 
nucleon). The quantity F in ( |13D is related to higher twist effects and it is suppressed by a 
factor A 2 /Q 2 . 

So, at high Q 2 and fi 2 » A 2 we get: 

F = CD, (17) 

where 

C = A(I-Vr( q ,k)\ k± ^ (18) 

with V and A defined above (see (0), (U)). 

Let us note that the quantity D (|16D is of the sum of the terms wich looks like the 
following (after an integration in transverse components of momenta fcj): 



z n—l k 2 



f d -^V{z n . 2 k ^)...( kl dkl f ^V(z ,^)D(z ,k 2 ,p 2 l (19) 

J Z n _2 "Vi— 2 J J Zq Ki 

where z = k n+ in/k n n. On the other hand, from an analogous formula in || we get the 
terms of the type 



J J Z n — l K n 

dz n _2 TT , kn-i\ f dzo TTt kf 



dk n _ 2 



Zn-2 



■U(z n ^^)...[^U(z ,^ 2 ), (20) 



where 



f/ = P At P /l V[/-(/-P /l )V]- 1 . (21) 

As one can easily see from (|T6|) , (|l9|) , P has the meaning of the distribution of partons 
whose virtualities vary ut to /i 2 . That is why we can consider C in Eq. fll7|) to be a coefficient 
function. It can be represented in another (equivalent) form: 

C = A[I-(I- K)V]-\q, k)\ k±=k _ Q . (22) 



Here 



K = P^P h . (23) 

So, the operator (I — K) — (I — Ph) + (/ — P^)Ph in Eq. (^) projects onto states without 
any collinear singularities. It acts on the full expression on the right [|TC|] : 

(/-(/- K)V)- 1 = / + (/- K)V + (/ - K)(V(I - K)V) + .... (24) 

The expression for Ci 1 ', calculated with the use of formula (|22"1), coincides with a s order 
result from Ref . |15| . 



Starting from Eq. (|18D or (^), we can represent DIS coefficient function C in the form 
analogous to (f|): 

C= (CO" -AC) (,,»)! (25) 



where 



= A(/ - V)- 1 (26) 

is a "naive" expression for virtual photon-parton coefficient function that does not take into 
account the evolution included in parton distributions D (0). In particular, if we consider 
a parton to be a gluon, we have C PhP = C PGF . 

From Eqs. fl25|), ( 22"D and ( |26|) we get the following expression for AC: 



AC = CKV(I ~V)-\ (27) 

where C(g, k) is given by formula (p2] ) but without imposing a condition k± = k^ = 0. In 
deriving (pTD we used a matrix identity 

(/ - y)- 1 -[/-(/- K)^}- 1 = [/-(/- K)l/] _1 iCy(J - I/)" 1 . (28) 

Let us note that both C PhP (q, k) (^B) and AC(q, fc)(P7|) contain, in general, singularities 
at A; 2 = 0, while their difference, C(g, fc) (^), does not. 



The formula (p7[) enables us to calculate AC 9 in any fixed order in a s . Let the coefficient 
function C and the kernel V have, respectively, the expansions: 

C = C (0) + C {1) + ... (29) 

and 

V = V (0) + V {1) + . . . . (30) 



Then we get: 



AC W = &?KVV> (31) 



and 



+ ^^« + C'( 1 )KV;«. (32) 



The kernels V can be related with parton splitting functions (a, b = q, g) flT4f : 

^-w^ = ^r 1 ' (^) , (33) 

where off-shell splitting function has the form: 



In particular, in the leading logarithmic approximation, we have: 

v£\ z y,e = o) = ^p!°\ Z ), (35) 

where pf£\z) is a leading order Altarelli-Parisi splitting function. 

If we put /i 2 ~ Q 2 in Eq. (^) and save (a s In Q 2 ) 2 and a 2 In Q 2 contributions, all the terms 
in expression (Qj for AC 9 can be reproduced (taking into account slight difference between 
our definitions of and AC 1 '"' and those from Ref. 0). Indeed, from ( p3|) we conclude 
that CPKVW ~ C« ®a s In Q 2 P§\ etc., where coefficient functions = C^ n \a s (Q 2 )) 
have no logarithms of Q 2 . 
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However, exact formula ( p^) results in additional contributions which are absent in (^j). 
In particular, due to power corrections in V^p (see (0)) the term CjP^ K[VjpVjp + V^Vjp] 
contains a non-leading contribution (a 2 InQ 2 ) in addition to a leading one ((a s In Q 2 ) 2 ). 

Let us note that expression ( p7|) and, consequently, (pTC]) and fl52|) do not have factorize 
forms as right hand sides of Eqs. (Q), do. Namely, the integration of C(q, I) or C( n \q, I) 
in the momenta I with the right part of the corresponding expression should be done. In 
particular, due to this integration in equation (|3lD we get a constant term in AC^ (g, /c) 
in addition to large logarithmic term. It depends both on a ratio f^ 2 /Q 2 and fi^/Q 2 , where 
/2q = — k 2 (simultaneously, we have to put k 2 = — /Xq in C PhP (q, k)). At /i 2 ~ Q 2 ^> m 2 and 
A; 2 = we get the constant term which is different from z(l — z) in (^). 

Fortunately, we are not forced to deal with the quantity AC (^7|) as we have derived 
the formula (p2|) which enables us to calculate C g and C c in any order in strong coupling 
without getting double counting in the coefficient functions. 

In the present paper we have studied the structure function F 2 . However, our formulas 
(P7|) may be also applied to the longitudinal deep inelastic structure function F L . 

In conclusion let us note that the problem of double counting should also exist for 
Wilson coefficients for light quarks, C q , considered in higher orders in a s . In such a case the 
double counting means that one and the same term is accounted for both in the coefficient 
function C q and in a distribution function D q . Formula (|18D ( (|22|) ) enables one to separate 
diagrams, describing C q , from those included in the evolution of the quark distribution D q . 
For instance, in the lowest order in a s we should get a term analogous to (f|). 

I would like to thank the Brazilian govemental agency FAPESP for financial support. 
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